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ABSTRACT 



Finite-amplitude standing wave effects in air at ambient 
conditions contained in rigid-walled cylindrical tubes with large 
length to diameter ratios were experimentally investigated. These 
results were compared to a perturbation solution of the one-dimensional 
non-linear acoustic wave equation which incorporates the dissipative 
effects due to viscous and thermal energy losses at the walls * The 
lowest resonance frequencies of the tubes ranged from 62.5 Hz to 
1 kHz, and sound pressure levels (based on the fundamental) ranged 
up to 160 dB. The finite-amplitude distortion was in excellent 
agreement almost up to the onset of shock. A detailed investigation 
of small amplitude attenuation in standing wave tubes was conducted 
and compared with the Kirchhoff equations. Agreement within one 
percent was obtained when consideration was given to the numerical 
analysis correction and end-effect losses. 
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1. Introduction 



It has been known for a long time that for the case of a plane 
elastic wave of any amplitude in a nondissipative medium the exact 
solution to the equations of hydrodynamics predicts a change in form 
as the wave travels ( 1 , 2 , 3 , 4 , 5 , 6) . In going one step further and 
considering that the medium is dissipative, it can be shown that only 
waves of relatively large amplitude evidence any appreciable change 
in form (4,7). This provides a good definition for finite-amplitude 
waves. Stated another way, a wave would be of finite-amplitude when 
the inclusion of second-order terms are required for the theory to 
predict a change in shape as the wave propagates. 

An important contribution beyond the work done by the early 
investigators of the nineteenth century was by Fay (8). He investigated 
the change of type in a periodic finite-amplitude plane wave of infinite 
extent in a viscous medium. He formed a Fourier-series solution valid 
far from the source that was an attenuated saw-tooth wave. 

In general the more recent theoretical approaches to finite- 
amplitude effects have been of two types. Approximation methods using 
some assumption on the manner of harmonic generation and absorption in 
a medium have been used (9, 10,11, 12 ^ and approximate solutions to the 
nonlinear differential equations by the perturbation technique or other 
means have been investigated (8,13,14,15,16). 

An expression for the spatial attenuation of an plane saw-tooth 
wave finite in extent was developed by Rudnick (17,18) using a "first- 
order” theory. Since the wave was confined, consideration was given 
to the attenuation due to the walls by using the results of Kirchhoff 
and Helmholtz. 
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Goldberg (13,14) undertook a solution of the problem of plane 
waves in an unbounded medium, by a direct perturbation approach on 
the equations of motion. He obtained a number of second-order 
solutions and was able to provide conclusions about the relative 
importance of nonlinearity and dissipation. 

The perturbation analysis has been extended to the sixth order 
by Keck and Beyer (15) for periodic plane progressive waves of infinite 
extent in a viscous medium. The wave equation is the same one treated 
by Fay, but the solution is valid near the source. 

A numerical solution for approximately the same problem has been 
obtained by Cook (20). An exact solution in parametric form was 
obtained by Blackstock (21) for the problem of a finite-amplitude 
sound wave produced by arbitrary motion of a piston in a lossless 
tube . 

Most investigators, like the ones already discussed, have confined 
their studies to traveling waves. However, a few have given consid- 
eration to standing waves. 

An exact solution of the nonlinear wave equation for the case of 
finite-amplitude standing waves in a closed tube as well as for one 
having a vibrating piston at one end was obtained by Keller (22). 

Upon assigning a value of -1 to the adiabatic exponent^? the nonlinear 
wave equation becomes linear, and this fact was used by Keller in his 
solution. These results "blow-up" at resonance because the dissipative 
forces are neglected in the solution. 

Betchov (23) constructs a solution at resonance based on discontinuous 
linearized solutions and a secular equation. This solution is interesting 
in that for an invisid fluid the amplitude at resonance is found to be 
finite and to be determined by non-linear effects. Betchov also discusses 
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the modifying effect of wall friction on his solution and assumes 
that this effect is equivalent to a body force proportional to 
the velocity. 

Saenger and Hudson (24) present a simple analytical description 
of the behavior of periodic shocks at resonance using a two part 
solution of the nonlinear one "dimens iona 1 wave equation. A point 
of importance in their solution is that the amplitude of the 
oscillations remains finite only because of the action of shear 
viscosity and heat conduction. 

Chester (25) has predicted asymmetries in finite amplitude 
standing waves with an elaborate theoretical investigation based 
on the nonlinear acoustic equations. The solution includes the 
effects of compressive viscosity, and of shear viscosity in the 
boundary layer at the walls of the tube; however, his solutions 
are not amenable to detailed experimental comparison. 

Coppens and Sanders (26) developed an extension of Keck-Beyer 
perturbation approach that included both the bulk losses and the 
wall losses. Their results provide information on the amplitudes 
and phases of the Fourier components of the distorted pressure 
waveform. The theory is applied to the problem of finite-amplitude 
standing waves in a cylindrical resonance tube that is excited by 
a piston at one end. This perturbation approach was extended by 
Ruff (27), and computer programs were written to extract more 
information from the theory. 

The purpose of this research was to conduct a detailed 
experimental investigation of finite-amplitude standing waves in 
air contained in rigid-walled tubes to verify the soundness of the 
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Coppens-Sander s theory as extended by Ruff. In addition, a detailed 
investigation of the attenuation of sound in rigid-walled standing 
wave tubes with attention given to end-effect losses was conducted 
to check the validity of applying the Kirchhoff attenuation equations 
to standing waves. 
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2. Theory 



2.1 General 

A plane elastic wave traveling in a fluid exhibits a change 
in form as it travels. This can be seen from the nonlinear properties 
of the equations of state and motion which follow. This theoretical 
development is extracted from that of Coppens and Sanders and is 
presented here to aid in interpreting the experimental work. 

The adiabatic equation of state of a fluid 

Cfj (!) 

can be represented by a truncated Taylor's expansion in the form 

CA/pP(p.p c ) + (e>/z/p)C/>-pS ( 2 ) 



where 



p = acoustic pressure 
O = instantaneous density at a point 
^> 0 = equilibrium density 

B/A = parameter of nonlinearity. 

If the fluid is an ideal gas, then B/A = -1 where ~0 is the ratio 

of principal specific heats. The coefficient A can be expressed as 



A "pA 



(3) 



where 



c^ = infinitesimal amplitude free field phase velocity. 



The equation of continuity for one dimensional wave propagation 
can be written in Lagrangian coordinates as 



pf/^p =/3 



(4) 
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where ^ is the particle displacement from the initial position 
at a. 

The force equation can be written in a slightly generalized 
form as 



+ (5) 

when u is the particle velocity and the operator F_ generates 

<-0 L 

those body forces in Lagrangian form that are active in the system 
of interest. 

The combination of Eqs . 2-5 yields 










(6) 



which for an ideal gas becomes 

3ef ~ ^ 0 ^ ' 3x§ ~ . ( 7 ) 

For acoustic processes where ^ Eq. 7 can be approximated by 

□;i km ') " v (8> 

Differentiation of Eq . 8 with respect to time yields 






[o L ■+ <£{ ) ci - ft Qz£) 



(9) 



where 

ft- / + B/A) 

. ci 

d)A r f ftftj 

k "A 2 - w 

Equations 6-9 have been the basis of a number of theoretical 
investigations utilizing successive approximations or perturbation 
expansions to describe finite-amplitude effects for traveling waves 
in fluids (with or without viscous dissipation) (8,9,10,15,21,28). 
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Under the assumption of weak dissipative processes and with the 

substitution of the perturbation series V ^*3 c 9 

o/„ y 

becomes 



y ja 2^1)%#' - & ■„ 



(id 



where the sum over j is understood to encompass just those frequencies 
appearing in (J) with 



**§*/**} 

0* ' 0 



♦ ‘n % ' 



/ 



W >/ 
sV */ 



(11a) 

(lib) 



To predict the finite-amplitude effects for one-dimensional wave 
propagation in a fluid system, it is only necessary to obtain the 
linearized, absorptive wave equation for the system, substitute the 
resulting expression for into Eq. 11, and solve for the frequency 
components . 

2.2 Rigid-Walled Tubes 

Wes ton 1 s (29) approach to the propagation of sound waves in gases 
contained in wide tubes can be used in the formulation of D , For the 



u 

wide tube the sound energy is diffused evenly over the tube cross 
section, and the adiabatic motion approximates to that in an unbounded 
gas. There is, however, a narrow boundary layer at the wall, analogous 
to the Prandtl boundary layer in continuous flow, where viscous and 
heat conduction processes lead to an energy loss. The losses in the 
boundary layer at the wall are much greater than the bulk losses. 
Weston* s equation in Euler ian form can be written as 



al & + f •-/) =■ <o, 

where _ r ^ w 

£ ~(G/S)CZt* 1 ' Z 



(12) 

(13) 
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r: ~f) 1 = Kirchhoff 1 



s constant 



(14) 



S = (uniform) cross-sectional area 
G = perimeter of area S 

CO = (angular) frequency of the propagating waveform 

$ = kinematic shear viscosity of the fluid 

i) /= thermal diffusivity or thermometric conductivity 

i = Y-Tl 

The real part of u represents the mainstream particle velocity 
in the tube. For Eq. 12 to be valid where ^ is the boundary 

layer thickness given approximately by 5 ('$/oJ) 2 (29). Equation 12 
can be written as 



t}^ u + C S'/co ) <$xgL ~ S iC —0 

with the assumption of one -dimensional flow. 

Then 



(15) 



(16) 



07) 



where ^ represents the value of the dissipation parameter ^ for a 
frequency • 

The approximate wave equation is now given by 

r f/jZj.AyL- . 4 ' -4- a , = .42£-f r/ ‘2*.: 

? L' r \ U OJj vaf ' 3°-' dCL ■■ 

“ C k h - 4 ; ; - 

For a resonance tube of uniform cylindrical cross-section excited at 
a = 0 by a piston with the same cross sectional area driven with constant 
acceleration at an angular frequency CO and terminated at a = L by a 
plane rigid boundary, the appropriate boundary conditions can be written 
as 

\ 



(0,6) ~ A 0 / a V ; o? (o> u ~ O)) 



(18) 
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{ Ljt) ~ .A-'O. 






^>/c (19) 



The phase angle ^ will be specified later. From these boundary 
conditions Coppens and Sanders present a solution that describes 
the finite-amplitude waveform. Certain quantities developed by 
Coppens and Sanders are important in any experimental investigation 
and are therefore reproduced here. The propagation constant, k, and 
the small amplitude attenuation constant, OC , are given by 



J, =^/e * 



( 20 ) 



and 




( 21 ) 



Q,will be discussed in greater detail in Section 3.3. 

Various frequency relationships in terms of the infinitesimal- 
amplitude resonance frequency, , are given by 

£tO CO -~CQp ( 22 ) 



for the frequency differential from resonance and 



U>4 





(23) 



for the frequency at which each harmonic making up the pressure wave- 
form has decreased in value to 1/ its maximum value. The 

dissipation parameter for higher order, , is given by 




(24) 



where the integer n represents the harmonic. It is also convenient 
to introduce the peak Mach number of the first-order solution 

M = U/ /c < 25 > 
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3 . Exper mentis 1 Considerations 



Insight into the finite-ampli tude case can be gained by first 
relating the theory to the infini tesiroal-ampli tude ca Sc* with a view 
toward isolating the experimentally determinable parameters . From 
the expressions 



and 



U, ■■ (^/cdXsm^kl *• SI 




(26) 



(27) 



the pressure ac the rigid end of the tube can be found to be 

.'4 

sfi' ( h /l- ff>C c A /&)($! Sik^CX-L ^Slri'Aj.) Si# 0)t n (23) 



When the acceleration of the piston is written as 



A fa. 6) - A 0 smfeoi 4 &cf^ } ) 



(29) 



it can be seen that fch piston Leads the first-order pressure at the 
rigid end of the tube by 90° - 0^ o The first-order pressure amplitude 
at the rig: - <-:t.d when the attenuation is small can be written as 

H, * U<i a a y [ ( « i- f * 5-. « a o \ " , 



(30) 



and the phase angle. 0 1 becomes 






C Q-//7 ■ A 



(of /,) 



l31) 



At resonance En;s 0 30 and. 31 become 



n 



l k 0 /ct^lC 



and 



f A 



> / it 

= (y 



(32) 



(33) 
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where the integer m represents the normal mode most strongly excited by 
the input frequency. When finite-amplitude effects are taken into 
consideration, the pressure is modified by the presence of the higher 
harmonics and can be written as 




It is convenient to introduce the resonance parameter Q and the phase 

n 

parameter 9^ defined by 



a 



/n 




(cOS 




(35) 



and 



'' ^ 

[can ^ 



= (4 -4, - 






(36) 



The dependence on the small amplitude attenuation coefficient is best 
introduced at this point by writing the dissipation parameter <f# } in 
the form 



J'l ~ / 






(37) 



through the use of Eqs. 24 and 21. 

If and 9^ are rewritten in terms of founder the assumption that 
the bulk losses are less than the boundary layer losses, then the 
approximate forms are given by 

!/? J.s ^ (38) 



Qa, <* Z CQ ^ Q 



and 



*#» % - / 



/* t 



K _ 



W.C'o). 



(39) 



For a resonance tube with a diameter of 1.125 inches, a length of 6 
feet, and a fundamental resonance frequency of approximately 100 Hz, 
the bulk losses do not equal the boundary layer losses until the input 
frequency i>s of the order of 10 Hz. This would be in the range of 
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the 10 harmonic, and it is readily seen that the assumption that the 
boundary layer losses dominate is valid. 

Equations 26-39 can be readily used for a comparison with the 

theory. The factors /^y^and c^ are known for the system to be 

used, and the quantities and are known input variables. 

O^cWq, can be determined experimentally. The pressure amplitude 

P can be measured directly, and 0 can be determined from the 
nn J n 

frequency dependence of the pressure measurements. 

Examination of certain of these terms reveals some important 
properties. Q is maximum above *and the higher the subscript n 
the higher the maximizing frequency. With constant acceleration, the 
higher harmonics have their maximum amplitudes at successively higher 
frequencies with narrower peaking, and the frequency region for most 
favorable finite-amplitude excitation lies between and cd x ~ U>p(C>/c)- 

The angle 9^ describes the difference in phase between the piston 
acceleration at a = o and the pressure at a = L. 

3.1. System and Equipment Description 

A diagram of the system used in this experimental investigation 
is shown in Fig. 1. The resonance tubes contained air at ambient 
temperature and pressure. In order to insure as nearly as possible 
ideal boundary conditions and to reduce the effects of mechanical 
coupling with the exciter, the resonance tubes were fabricated of 
steel with an inside diameter of 2.250 inches and a wall thickness 
of 1.125 inches. The inside surface of the tube wall was cleaned, 
but no other attempt was made to improve the finish. The rigid end 
was terminated with an aluminum cap 2.0 inches thick. To insure no 
leaks, the surface of the cap and the end of the tube was machined 
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and polished and high vacuum silicone grease was used as a sealer. 

The cap was held in place by five studs screwed into the tube wall. 

A total of five tubes (identical except for length) were used to 
give a spread in fundamental resonant frequencies in excess of an 
order of magnitude. The lengths of the tubes were 7, 14, and 28 
inches and 6 and 9 feet. The six foot tube was used for the major 
portion of the investigation and for comparison with the computer 
predicted results. 

The piston design was critical for several reasons. It had to 
be heavy enough that a minimum of distortion was introduced by the 
force exerted upon it by the distorted standing wave, yet light enough 
that desired acceleration amplitudes could be obtained. A piston 
weight of approximately 0.8 pounds was used. A clearance of 
approximately 0.006 inch was provided between the seat and the piston, 
and an air-tight seal was obtained by the use of two 0-rings and 
silicone grease. The 0-ring grooves in the piston were designed so 
that the neoprene 0-ring would be under compression and were spaced 
to provide a reservoir for the grease. A smooth seat for the piston 
was provided by machining and grinding the inside of the tube. The 
plane surface of the piston was machined and polished. 

The M-B Electronics Model EA 1500 exciter due to its versitility, 
small size, and maximum acceleration rating of 124 g*s was chosen 
as the acoustic driver. For the production of intense standing waves 
in a resonance tube, large displacement amplitudes are not needed, 
but large acceleration amplitudes are needed. The permanent magnet 
vibration generator with its large force ratings and small moving 
mass readily meets these requirements. The frequency and amplitude 
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can be varied simply by adjusting the oscillator output. The EA 1500 
has a frequency range from 5 Hz to 10 kHz which is well beyond the 
first cross mode of operation and the limiting value for this investi- 
gation. 

The M-B Electronics model 2120MB power amplifier was selected to 
complement the exciter. The distortion was less than 0.77, at rated 
output of 125 watts in the desired frequency range of operation and 
less than 0.37, at driving levels below shock. A single unit provided 
accelerations to 50g*s for the piston used. Higher levels of acceleration 
can be obtained by paralleling additional power amplifiers. 

The severe requirements for a very highly stable oscillator with 
a precisely repeatable frequency setting led to the choice of the 
General Radio Type 1161A-A6C coherent decade frequency synthesizer as 
the primary frequency source. The synthesizer has a readily repeatable 
frequency setting to better than 0.001 Hz. The synthesizer was checked 
repeatedly with a counter, and no variation in the frequency output was 
measured. The output level control was modified to permit finer adjust- 
ment. A General Radio Type 1160-P2 Sweep and Marker generator was used 
in conjunction with the synthesizer. 

The piston motion was monitored by the use of a Endevco model 2215 
accelerometer mounted inside the piston. This model will work directly 
into a VTVM and still provide satisfactory frequency response up to 
10 kHz. It was necessary to mount the accelerometer on the centerline 
of the piston and directly to the leading component of the piston to 
reduce distortion and improve the accuracy of the results. 

The Hewlett-Packard Electronic Counter Model 521C was used to 
monitor the frequency input and the frequencies from the wave analyzer. 
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No deviations were noted between the set frequency and the output 
frequency for the synthesizer. 

A Bruel and Kjaer Type 4136 condenser microphone was mounted 
in the rigid end so that the diaphragm was flush with the end of 
the tube. T he small diameter (only % in . ) was chosen to minimize 
the disturbance of the rigid boundary condition. The frequency 
response of the microphone was flat to within +0.5 dB from 80 Hz 
to 20 kHz, well beyond the range of interest for all five tubes. 

The rise time of this microphone was 4.5 microseconds. 

All voltage measurements were taken with Hewlett-Packard 400D 
Vacuum tube voltmeters, and the waves shapes were observed using a 
Tektronic Type 565 dual-beam dual-trace oscilloscope. The distortion 
present in the wave forms was analyzed for harmonic content using 
a General Radio Type 1900-A wave analyzer. This analyzer has a band 
width selection of 3, 10, and 50 Hz. Some Q-curves were plotted 
using data Equipment Company Model 800 AM X-Y recorder and the 
General Radio sweep and marker generator. 

Figure 2 is a schematic of the system used to monitor the 
temperature of the gas. Because of its linear resistance character- 
istic in the range of interest, a Fenwal Electronics Inc. Model 
GA 51P1 thermistor was used and was mounted in a small hole drilled 
in the rigid end plate. The leads were insulated, and the unit was 
set in epoxy with the glass bead at the surface. The epoxy provided 
an excellent seal and a means of smoothing the end plate to reduce 
distortion. The temperature read-out was provided by means of a 
microammeter; when a continuous record was needed, a strip recorder 
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was used. The system was calibrated with the aid of temperature 
baths of acetone and dry ice, ice and water, and boiling water. 

With this system, temperature changes in the tube of less than 0.1 C° 
could be easily detected. 

3.2. Calibration of Equipment 

The experimental determination of P requires that the absolute 

nn 

calibrations of the microphone be known. In addition, the determination 

of 0 /j requires that the absolute calibration of the accelerometer also 

be known. The accelerometer sensitivity is defined as S A = V A /A where 

A A o 

is the peak amplitude of the accelerometer output voltage, and the 
microphone sensitivity is defined as = V^/P where V^. and P are the 
peak amplitudes of the microphone output voltage and the pressure 
respectively . 

The calibration of the accelerometer was accomplished by direct 
measurement of the displacement amplitude of the piston face for a 
range of accelerometer output voltages. Displacement amplitudes up 
to 0.060 cm were measured with a reproducibility of 0.0002 cm using 
a traveling microscope and a stroboscopic light. At 100 Hz, fifteen 
determinations of the accelerometer sensitivity resulted in 
£>A - (6.66+ 0.03) x 10 ^ vol t/ (m/sec^) . 

The microphone was calibrated using a Bruel and Kjaer Model 
4220 pistonphone which produces a known acoustic pressure level of 
124 + 0.2 dB at a frequency of 250 Hz. The sensitivity was found 
to be S M =(-76.4 + 0.2)dB re volt/^b = (1.50 + 0.03) x 10" 3 volt/(Nt/m 2 ) . 
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3.3. Determination of the Attenuation Constant 



The necessity of having an accurately determined small amplitude 
attenuation coefficient and the reported wide variation in agreement 
between the theoretical and the experimentally determined attenuation, 
particularly in tubes containing standing waves, prompted a more 
thorough investigation of this problem. 

3.3.1. Background 



the effects due to viscosity, the causes of the sound attenuation 
in cylindrical tubes have been widely discussed and investigated 
(29,31-42). Experiments of Kundt led to Kirchhoff's (6) complete 
theory for the wide tube (defined in Section 2.2) which considered 
both viscosity and heat conduction. In general, investigators 
report that the Kirchhoff equation gives the proper frequency, 
radius, and pressure dependence but that the experimental absorption 
is greater than calculated. Weston (29) presented more accurate 
expressions for the velocity and attenuation by including more terms 
in the expansion of the Bessel functions. Weston's expressions for 
velocity and attenuation for the wide tube (which we are considering) 
are : 

Phase velocity of sound: 



Beginning with Helmholtz in 1863 (30)^ who considered only 




(40) 




) 



cm -ssg 



Attenuation : 




P - 



r 'j 



(41) 
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where 



tj 7 - YT 1 



(42) 



/ 



W) O'/I)*- --0S74- at,/s«!\ 



The first term of Eq . 41 is the ordinary wall absorption, the second 
the ordinary free gas attenuation, and the third arises because the 
sound energy is not uniformly distributed over the tube. The 
numerical values here are for dry air a 20°C and 760 mmHg and were 
checked using constants obtained from National Bureau of Standards 
Circular 564 (43). 

Shields, Lee and Wiley (44) present a numerical solution of the 
Kirchhoff equation and state that Kirchhoff s phase velocity correction 
is always too large (giving values which are too small), and that the 
calculated absorption is always too small. They also give a method 
for determining the error in the Helmholtz-Kirchhof f equations. 

The work thus far presented is for progressive waves in infinitely 
long tubes. It has been stated by Shields (44) that neither the 
Kirchhoff tube absorption and velocity equations nor the numerical 
analysis would be expected to apply to standing wave tubes. 

In considering resonance tubes, there is almost universal agreement 
on the existence of excess experimental absorption (44-48) when compared 
with the absorption calculated using the equations for traveling waves. 
For example, Parker (45) reported an excess which varied from 6.5 % 
for helium to 30% for argon; Smith and Tempest (44) report 30%; 

Henderson and Donnelly (46) found 107o. 

3.3.2. Theory for Standing Waves 

For standing waves to be generated in a rigid-walled tube, 
an integral number of half wave lengths must be contained in the tube 
length, i.e., ^ ^PL//\y\ where m = 1,2,3,... From the fact that 
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yt = c/f 
ys m m 



/?) 



'>Ti 






(43) 



where c is the velocity of sound in the fluid contained within the 
tube. For a non-dissipative medium, the infinitesimal pressure at 
the rigid end (where a = L) is given by 




OOS 



£/~J (:> 



) / (so s'^o - j&L ) 



* 



(44) 



If the medium is dissipative, the sound pressure at the rigid 
end is still given by Eq. 44, but now k may be considered as complex 
in order to account for the dissipation. We can now write Eq. 44 as 



i > 

r 



(/o^/4 0 /co)[smtr&L cos^t (^) 

which is the same as Eq. 28 with a 90° phase change in the acceleration. 
And if the absorption is small (i.e., L<*< i ) , the amplitude of Eq. 45 

is given by Eq. 30 

,i . v r- . a ^ -K 

(30) 



-T'/l 
c 



p - (fip 0 Aa/cp)[C(/.LT 1 

This is of the form of a typical resonance curve if the width of the 
resonance curve is much less than the fundamental frequency. From this 
resonance curve, the half-power, or the "down 3dB" points, can be used 



to calculate the absorption in the resonance tube by 






C 



(46) 



This method for determining the attenuation was used by Parker and 
by Henderson and Donnelley. 

For the resonance tube, end effects seem to be the major additional 
consideration. If the losses at the ends are considered to be due only 
to heat conduction and are not a function of the viscosity, then an 
expression for these losses may be obtained by using the heat conduction 
portion of li* . 
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The relationship of the viscous and heat conduction terms is given by 



Eq. 42 



*j/ 






2> " ’ ' ’ ' " ' ' ' ' (42) 

where the first term is the viscous term and the second term is the 
heat conduction term. Assume that the attenuations are additive, 
i .e . , that 

rl 1 - cC 6 -/-o: e (47) 

whereC^ is given by Eq. 41 and (X. is due to the losses at the end. 
Kirchhoff's expression for the attenuation due to the walls alone is 

°V ' O', cvO 



which becomes 



a- 



= ^ 97V 'r 

Oj y ' < (A 



(48) 



(49) 



when evaluated for dry air at 23 C and atmospheric pressure. Equation 
49 for tubes with a radius of 1.125 inches becomes 



c(* M - x /o '\£ ,/2 



(50) 



Assuming that the end has only heat conduction losses, then 

cC e * (51) 

Since the heat conduction term is approximately one-half the viscous 
term for dry air at 23°C and atmospheric pressure, the attenuation due 
to the ends can be related to the attenuation due to the walls by 

O^e ^ C(. OJ ^i%-/dL (52) 



3.3.3 Experimental 

Attenuation constants were experimentally measured for five 
identical resonance tubes of varying length by two methods. The first 
consisted of measuring the width of the infinitesimal-amplitude 
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resonance curve between the half-power points and using Eq. 46 to 
calculate the attenuation constant. This will be referred to as 
the Q method. The second method consisted of measuring the 
microphone and accelerometer output voltages at the resonant 
frequency and using Eq. 32 which can be written in terms of 
voltages 



1/ W /oZ-exp rrr ^ I Sy, / $.%) (53) 

A plot of vs V^, for , should result in a straight line 

of slope ( A /(?' /^ 7 ^D(S /S ). From this slope, (% can be calculated 
• ' ii A exp 






if S.. and S A are known. This method will be referred to as the A/P 
M A 



method. 



3.3.4. Comparison of results 

Attenuation constants were calculated for all five tubes using 
Eq. 49 and Shield’s corrections for dry air at 23°C and 760 mmHg. 

These values for CX co j are compared with the experimentally determined 
(j'/C for the fundamental modes in Table 1. The experimental results 
for all five tubes by both methods are compared with the theory in the 
range of interest in Figs. 3 and 4 where the log of the attenuation is 
plotted vs the log of the frequency. There is almost perfect agreement 
between the calculated and the experimentally observed attenuation 
constants by the Q method. This is true through the 50th harmonic 
for the 9* tube and the 40th harmonic for the 6' tube. At these 
points side modes were calculated to be present and were observed by 
distortion in the Q-curves. 

Except for the 7 in. tube, the effect of the end- loss correction 
is too small to see on the graphs. As a check on the end-loss 
approximation, if + <X Q and OC in Table 1 are compared, 

it is seen that the assumption that the end losses are due solely 
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to heat conduction is slightly low. If the same losses are 

assumed for the ends as for the side wall, then an R/L relationship 

would hold. Again compare CC m and rY in Table 1, and it can 

T exp 

be seen that this correction is too large indicating that the 
end correction is somewhere between these two extremes. If an 
empirical relationship of R/2L is taken, then the agreement is 
better than 1% for all but the seven inch tube. 

The results of these investigations indicated that the Kirchhoff 
equation (using Shield 1 s corrections and properly corrected for end 
losses) can be used to determine the attenuation in standing wave 
tubes of high Q where the boundaries approach ideal boundaries and 
where the length of the tube is much greater than the radius (on 
the order of ten to one) . 

The results further indicate that both the Q method and the A/P 
method are valid methods for determining the attenuation. The Q 
method has the advantage of requiring very little data and data 
processing to obtain the attenuation when a highly accurate frequency 
source is available. However, the acceleration had to be closely 
controlled, or an attenuation constant greater than actually existed 
was found. This change in the impedance of the driver near resonance 
could account for some of the problems in the past. The A/P method 
is reliable and useful when the frequency is not critical though the 
accelerometer and microphone must be accurately calibrated. 

The effects of temperature cannot be minimized. A 3°C change 
makes approximately a 1% change in Kirchhoff* s constant for dry air 
as given by Eq. 42 and a resulting change in (X . Also of interest 
is that a 0.6°C temperature change will result in a 0.1 Hz change 
in cOp for the six foot tube. Therefore, it was important that the 
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temperature of the system be closely monitored, particularly when 
the resonance tube was being excited at very high levels. 

3.4. Determination of the Resonance Frequency 

The resonance frequency was determined by finding the frequency 
at which P- was a maximum for a fixed A . Since ^O yr . was defined on 
the basis of the first-order solution, was set low enough that 
finite-amplitude effects were negligible. Repeated measurements 
indicated that Cty* could be determined to better than 0.01 Hz 
providing consideration was given to temperature. As a check on oCfa , 
the frequencies of the ,f down 3 dB" points (which could be determined 
to within 0.001 Hz) were determined, and the mid- frequency was compared 
with Oj> ¥ . The agreement was within 0.005 Hz. 

To arrive at a value of characteristic of the conditions 

actually prevailing at the time that the finite-amplitude data were 
collected, the experimental procedure was designed so that 10, n and the 
temperature were determined before the run and after the run and the 
temperature was recorded throughout the run. This technique insured 
that ^ was known to within approximately 0.02 Hz for all runs. 

3.5 System Alignment 

It was necessary that the piston be very carefully aligned before 
each run or distortion was introduced into the input waveform. The 
system was designed so that the resonance tube was fixed and the shaker 
was mounted on an vibrational insulated base plate that was provided 
with three dimensional adjustment. 

The alignment procedure consisted essentially of adjusting the 
piston in the tube until it was free to move . The piston was then 
excited ata low level, and the accelerometer output was monitored on 
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one trace of a dual trace oscilloscope. A sine wave from the 
frequency synthesizer was displayed on the other trace. Fine 
adjustment was made to the shaker base until the two wave forms 
appeared to be identical. The excitation level was then increased 
to the maximum level desired and additional fine adjustment was 
made as necessary. 

The distortion present was determined by measuring with a 
wave analyzer the percentage of harmonic content relative to the 
fundamental in the accelerometer output. Due to the impedance 
mismatch between the accelerometer and the wave analyzer, it was 
necessary to use some form of impedance matching, and a Burr-Brown 
Model 100 preamplifier was used because it had a low noise and 
provided excellent isolation characteristics. 

By using the visual method of waveform comparison for alignment, 
it was found that the second harmonic distortion of the acceleration 
waveform could be kept below 0.27o of the fundamental and the third 
harmonic distortion could be kept below 0.37o of the fundamental. 

All other harmonics were masked by the noise level. 
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4. Measurements and Analysis 

Three measurement procedures were employed to compare the 
experimental results with the theory. These procedures provided 
information to compare the pressure waveforms, Q-curves, and a 
Fourier analysis of the waveforms with these quantities predicted 
by theory. 

4.1. Pressure Waveform 

The pressure waveforms for various sound pressure levels were 
compared with those predicted by the computer program FINAMPI. The 
program FINAMPI computes and graphs the predicted pressure waveform 
at the rigid end using all terms in the perturbation solution through 
the sixth order with input variables of ay and M. 

The measurement procedure consisted of first setting A^ and the 
frequency at either a ) or at some desired frequency off resonance. 

The pressure waveform was then monitored on one channel of a dual- 
beam, dual-trace oscilloscope. A sine wave from the frequency 
synthesizer was displayed on the other channel for reference and 
distortion determination. Photographs were made of the oscilloscope 
trace for comparison purposes. 

The analysis consisted simply of visually comparing the photographed 
pressure waveforms at the various SPL’s and frequencies with those 
predicted by the computer program. 

4.2. Fourier Analysis 

A Fourier Analysis of the pressure waveform was made. The same 
procedure was used here as with the waveform analysis. In addition the 
harmonic content was measured as for the Q-curves. 
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The analysis consisted of applying a grid to the photograph 
and taking the amplitude of 64 equally spaced points in one cycle* 

These data points were used as the input for the computer program 
FOUANAL . FOUANAL takes these data points and does a Fourier snythesis 
on the wave shape* The amplitudes of the Fourier components were 
then compared with the measured values of the harmonics. 

4.3. Q-Curves 

It is possible to relate the response of the resonance tube to 
the frequency response of an electrical network. In this case the 
pressure curves of the resonance tube would be analogous to the "Q" 
response curves of the electrical network. Ruff has developed an 
expression for the amplitude of each harmonic of the pressure waveform 
as a function of ZVx> and M. The procedures outlined here provided 
the means for displaying the harmonic content of the pressure wave 
as a function of frequency near resonance for a desired value of A q . 

The measurement procedure consisted of first setting the value 

of A and then measuring with a wave analyzer the pressure of the nth 
o 

harmonic, P^, at the rigid end as a function of frequency. The 
frequency was swept from well below the half-power point to well above 
the half-power point for the fundamental in steps of 0.02 Hz. For the 
other harmonics through the sixth, the frequency was swept from just 
above the noise level through the area of interest back to just above 
the noise level in steps of 0.01 Hz. 

The data were analyzed by first plotting the various P in term 
of voltage versus the frequency. The data were then normalized using 
the magnitude of the leading term of the particular harmonic taken 
from the computer program QCURVE and QCURC. This was done to facilitate 
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comparison with the normalized graphs that were drawn by the computer. 
The fundamental resonance frequency was determined and used as the 
reference frequency for comparison with the computer graphs. The 
normalized data points for each harmonic were plotted and compared 

to the computer graphs provided by the two programs, QCURVE and QCURC. 

QCURVES is a program that uses the perturbation solution and is 
designed to compute and graph the resonance response of each harmonic 
using terms of the sixth order as the fundamental is swept through 
resonance. QCURC is a program similar to QCURVES except that it uses 
the leading terms of each harmonic in the perturbation solution. 

The major points considered were the amplitude of the harmonic, 
the frequency at which the maximum amplitude for each harmonic occurred, 
and the bandwidth of the various harmonics. 
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5. Results 



In this section the experimental results will be compared with 
the theoretical results obtained by Ruff. 

5.1. Pressure Waveforms 

The first comparisons made will be of the pressure waveforms 
(Figs. 5-8). In these figures the large drawing is of the waveform 
predicted by the perturbation approach (program FINAMPI) and the 
smaller inset is a drawing from an oscilloscope photograph of the 
actual pressure waveform. The dashed line in both drawings is a 
sine wave which is included as a reference. 

In the low Mach number region which includes Figs. 5 and 6, it 
appears that the theory and experiment are in very close agreement. 

In this region there is relatively little finite-amplitude distortion, 
and good agreement is expected. When the Mach number exceeds a value 
of about M = 0.005, the finite-amplitude distortion reaches a 
significant level. However, Fig. 7, which is for a Mach number of 
M = 0.007 indicates that the theoretical and experimental results 
are still in surprising agreement. As the Mach number is increased 
into the realm between significant finite- amplitude effects and the 
region of shock, the theory begins to breakdown as is evidenced in 
Fig. 8. This figure is for M = 0.009, and it is seen that the 
theoretical waveform displays details which are not observable in the 
experimental waveform. 

It is evident that the theory becomes invalid when it is extended 
into the region of relatively high Mach number. 
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5.2. Fourier Analysis 

In addition to the foregoing, the experimentally determined 
values of the amplitude of each harmonic may be compared with the 
results from the program FOUANAL. 

The results of the comparison with the program FOUANAL are 
included in Table 2. FOUANAL uses the data taken point by point 
from an oscilloscope photograph and computes the Fourier coefficients 
of the various terms of a Fourier expansion of the waveform. The 
small size of the photographs and the thickness of the trace line 
in the photographs were such that only qualitative results were 
expected. 

FOUANAL also computes the phase angle associated with each 
harmonic. These results are included in Table 3. 

5.3. Q-Curves 

Q-curves for three Mach numbers (M = 0.004, M = 0.005, M = 0.009) 
were experimentally obtained and graphed for comparison with the Q- 
curve program results for the same Mach numbers. These computer program 
results indicated by a dashed line for QCURVES and a solid line for 
QCURC are contained in Figs. 9 through 26. Several things should be 
mentioned concerning these graphs. First, the three Mach numbers chosen 
cover the range from very slight distortion to significant distortion 
and should indicate the usable extent of the theory. Second, in every 
figure, the Q-curve results from QCURVES, lies below the results 
obtained using QCURC (the first term method) . The data points in 
general fall very close to the results of QCURC except for very high 
Mach numbers where the predicted amplitude of the higher harmonics is 
much larger than the experimentally determined. Another interesting 
fact is that the same harmonics for each Mach number have exactly the 
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same Q-curve using the first term method. This is not surprising 
when it is remembered that in the normalization scheme the dependence 
of the leading terms will not change for different Mach numbers. 

5.4. General 

It would have been interesting to obtain the same results for 
the other tubes used in the attenuation measurements , but time did 
not permit. However, some observations that were made are worthy of 
note. The distortion of the pressure waveform developed in the same 
way for all five tubes; however, in the formation of shock slight 
differences were noted and were more prominent in the short tube. 

A sequence of photographs of the pressure waveforms for the 
14 in., 6 ft., and 9 ft. tubes at sound pressure levels of 155, 157, 
and 160 dB were taken and are shown in Figs. 27, 28, and 29. 

The wave shapes are very similar; however, there are certain 
noticable differences. The small spur at the base of the shock is 
present in all waveforms, but is more prominent in those for the 6 ft. 
tube. This point needs further investigation to determine whether 
the pressure actually decreases immediately before the shock or if 
the spur is due to some property of the instrumentation. The steepness 
of the shock front increases as the resonant frequency increases. The 
grass in the pressure waveform increases with an increase in resonant 
frequency and the drive level. This is a result of the "ringing 11 in 
the piston introduced by the shockwave banging against the piston face. 
The "grass" is also more prominently displayed at the higher frequencies 
due to the relative time scales and the higher drive levels required to 
attain the same sound pressure level. At very high drive levels, the 
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shock wave striking the piston can be observed in the accelerometer 
output from the piston, and its position can be shifted by small 
changes in the drive frequency a 



6. Conclusions 



Experimental investigation has shown that the theory accurately 
predicts the position relative to the fundamental at which the maximum 
amplitude of each harmonic occurs throughout the range from very low 
to very high Mach numbers. The investigation has also shown that in 
the region of relatively low Mach number the theory accurately predicts 
the amplitude of the first six harmonics that comprise the pressure 
waveform. The fact that the third harmonic is slightly high at the 
lower Mach numbers can be accounted for in part by the third harmonic 
distortion in the input. At high Mach numbers the predicted amplitudes 
of the higher harmonics are much greater than the experimentally 
observed ones. At this writing there is no definite explanation of 
this . 

The reconstructed pressure waveforms drawn by the computer are in 
very close agreement with the experimentally observed ones at relatively 
low Mach numbers. When the finite-amplitude effects become pronounced, 
the theory does not properly predict the spur that forms just prior to 
the shock formation; however, the theory does indicate a definite 
distortion of the waveform in this region. 

The Fourier analysis of the pressure waveform shows good agreement 
with the experimentally determined amplitudes at all Mach numbers. 

The experimental phase angles show an excessive scatter, and the 
agreement with the theoretical values is* only qualitative. This appears 
to be the only place where the experimental results are less reliable 
than the theory. Since the major part of the distortion takes place 
in a relatively small portion of the waveform, the accuracy of the 



44 



Fourier analysis could be improved through the use of more data points 
in the region of the developing shock. It is suggested that further 
work along this line be postponed until the oscilloscope waveform can 
be more closely linked with the computer input. 

In addition, the dependence of the theory on the infinitesimal- 
amplitude attenuation coefficient has been demonstrated. 

In general, then the statement that the perturbation theory can 
adequately describe the effects of finite-amplitude standing waves 
in rigid-walled tubes has been shown to be valid provided the Mach 
number is such that the distortion is not too great. This of course 
corresponds to the pre- shock area, and even then precludes much of 
that area in which significant distortion is evident. 
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2.45 


-1.6 


7" 


3.255 
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a 



a 
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O' 



O' 



exp 



Attenuation constant using Weston's formula 
<y + Shield's correction 

Attenuation constant due to the end effect 
%al* r/ 3L or "cal* r/ L 

Experimentally determined attenuation 



ATTENUATION CONSTANTS AND END-EFFECT CALCULATIONS FOR 
FIVE RIGID -WALLED TUBES CONTAINING STANDING WAVES 

TABLE 1 
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APPARATUS AND INSTRUMENTATION 





Temperature Monitoring System 
Figure 2 
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Pressure Waveforms 
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Figure 




Q-curves of the Fundamental Harmonic 
M = 0.004, SPL = 147 dB, L = 6' 



Figure 9 
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M = 0.004, SPL = 147 dB, L = 6 1 



Figure 10 
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Q-curves of the Third Harmonic 
M = 0.004, SPL = 147 dB, L = 6* 



Figure 11 
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Q-curves of the Fourth Harmonic 
M = 0.004, SPL = 147 dB, L = 6' 

Figure 12 
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Q-curves of the Fifth Harmonic 
M = 0.004, SPL = 147 dB, L = 6' 

Figure 13 
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Q-curves of the Sixth Harmonic 
M = 0.004, SPL = 147 dB , L = 6' 

Figure 14 
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Q-curves of the Fundamental Harmonic, 
M = 0.005, SPL = 149 dB, L = 6* 

Figure 15 
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Q-curves of the Second Harmonic, 
M = 0.005, SPL = 149 dB, L = 6* 

Figure 16 
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Q-curves of the Third Harmonic, 
M = 0.005, SPL = 149 dB, L = 6' 

Figure 17 
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Q-curves of the Fourth Harmonic, 
M = 0.005, SPL = 149 dB, L = 6* 

Figure 18 
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Q-curves of the Fifth Harmonic, 
M = 0.005, SPL = 149 dB , L = 6 1 

Figure 19 
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Q-curves of the Sixth Harmonic, 
M = 0.005, SPL = 149 dB, L = 6' 

' Figure 20 
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Q-curves of the Fundamental Harmonic, 
M = 0.009, SPL = 154 dB, L = 6' 

Figure 21 
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Q-curves of the Second Harmonic, 
M (= 0.009, SPL = 154 dB , L = 6* 

Figure 22 
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Figure 23 
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Q-curves of the Fifth Harmonic, 
M = 0.009, SPL = 154 dB, L = 6' 

Figure 25 
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Figure 26 
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Pressure Waveforms 
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